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1 Introduction 

The study of two-dimensional lattice models in Statistical Mechanics advanced dramat- 
ically with the advent of Baxter's concept of commuting transfer matrices [l]. This 
method introduced the concept of integrability in Statistical Mechanics and paved the 
way for the development of a variety of exact methods exploring the aforementioned 
commutativity. As examples of those methods we have Baxter's T — Q relation [l], the 
algebraic Bethe ansatz |2||, the inversion trick |3], etc. Also a variety of models can be 
tackled by those same techniques such as vertex models, solid-on-solid models and hard 
square models. 

Nevertheless, the implementation of those methods depends drastically on the bound- 
ary conditions chosen and the case of domain wall boundaries deserves special attention. 
This kind of boundary condition for the six vertex model was introduced by Korepin in [2] 
who also obtained a recurrence relation determining the model partition function. This 
recurrence relation was later on solved by Izergin in terms of a determinant fs]. More- 
over, the exact solution of this model raised the issue of the sensitivity of the six vertex 
model bulk properties with respect to the boundary conditions in the thermodynamical 
limit [6]. 
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Figure 1: Two-dimensional lattice with M x N retangular cells. 

A natural question that emerges in this scenario is how this sensitivity with respect 
to boundary conditions extends to the eight vertex model. In that case, although the 
partition function has been evaluated in (^-[9], the lack of manageable expressions have 
eluded the analysis of the thermodynamical limit except for a particular value of the 
anisotropy parameter 
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Motivated by this scenario and keeping in mind the relation between Baxter's eight 
vertex model and solid-on-solid models, also refereed to as SOS models, here we demon- 



strate that the algebraic-functional method introduced in 11 , 12 can also be used in that 



case. Using that method we derive a multiple integral formula for the partition function 
of the trigonometric SOS model with domain wall boundaries. 

This paper is organised as follows. In the Sec. [2] we briefly describe SOS models in 
Statistical Mechanics with emphasis on the case of domain wall boundary conditions. In 
the Sec. [3] we present the dynamical Yang-Baxter algebra and demonstrate how it can be 
explored in order to obtain a functional equation for the partition function of the SOS 
model with domain wall boundaries. This functional equation is analised in Sec. |4] and 
the solution of the equation is given in Sec. |5] as a multiple contour integral. Technical 
details are presented in App. [K\ through App. |F] 



2 Solid-on-solid (SOS) models 

We consider a two-dimesional lattice formed by retangular cells juxtaposed as in Fig. [T] 
For a lattice with + 1 rows and M + 1 columns we have M x N retangular cells 



and we associate the Boltzmann weight wij ( "^1' '■' ''*+^'-?+^ J to the cell enclosed by the 

cartesian coordinates + 1), {i + 1, j) and (z + 1, j + 1). Each retangular cell 

is simply referred as face and its configuration is characterised by the set of variables 
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Figure 2: Face at coordinate and its Boltzmann weight. 

+ h + 1,2 ... ... ... li + l,M+l 



Figure 3: Face monodromy matrix. 



{h,j, k,j+i, U+i,ji ^i+ij+i}- Diagrammatically this association is depicted in Fig. [2] and the 
partition function of the system is then given by 



N M 



li 



k+i,j 



i+1 
+1 



(2.1^ 



Face monodromy matrix. Let us define a matrix T with components 



/T-'i + l,lvi'i + l,Af + l 
k,l,---,k,M+l 



M 

n 



Wi 



h+i,j h+i,j+i 



(2.2) 



We shall refer to T as face monodromy matrix and it is diagrammatically represented 
in the Fig. |3| For simplicity of notation we also refer to the components of the face 
monodromy matrix defined by ( |2.2[ ) as 7^*^^. Employing this notation the partition 
function Z reads 

Z=J2 ^"'^^ ■ ■ ■ '^N^' ' (2-3) 



{k,j} 



and the lattice in Fig. [l] can be entirely built as the product of elements Tl^"" representing 
each horizontal layer. 



Boundary conditions and integrability. Following Baxter 13 we consider the face 



monodromy matrix (2.2) with periodic boundary conditions in the horizontal direction. 



This corresponds to setting U^m+i = and under these conditions it is convenient to 
define an operator T with components 



t; 



•k+l,l,---,k + \,M 



k,i^---ik,M 



, h+1,1 h+l,2\ f h+1,2 h+1,3 



li+l,M h+l,l\ ^2.4) 
H,M HA J 
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Figure 4: Graphical representation of the Yang-Baxter relation for SOS models. 



The operator T is commonly denominated (face) transfer matrix and it plays an im- 
portant role in estabilishing the integrability of two-dimensional lattice models p3] . We 



proceed by defining a second transfer matrix T' similarly to (2.4) but with Boltzmann 
weights w' instead of w. Next we look for conditions on w and w' such that the transfer 
matrices T and T' form a commutative family, i.e. [T, T'] = 0. This requirement leads 
to the following relation 



E 



h lo\ fh h , 



^6 ^5 
lo I4 



E 



h lo 



lo h 

1 7 I '^u+v I J I I '^v \ 1 1 



h k 



(2.5) 



where u and v are complex variables parameterising the manifold where the transfer 
matrices form a commutative family. For a detailed derivation of (2.5) we refer to [M]. 



The Eq. (2.5) is usually referred to as Yang-Baxter relation 15 , or simply Hexagon 
identity, and the local equivalence transformation described by (2.5 ) is depicted in Fig. |4 



Moreover, if we also consider periodic boundary conditions in the vertical direction, i.e. 
In+ij = li,j, the partition function (2.1) becomes simply 



Tr (T^) 



(2.6) 



and its evaluation is translated into an eigenvalue problem 16,17. Throughout this 



paper we shall consider a different class of boundary conditions where (2.6) does not 



apply, although the bulk model is still governed by statistical weights satisfying (2.5). 



The trigonometric SOS model. The variables {/jj, /jj+i, /j+ij, Zj+ij+i} depicted in 
Fig. |2] are also called height functions and they characterise the configuration of the 
associated face. The degree of freedom lij is also refereed to as the colour of the face at 
the position Furthermore, we can also impose restrictions on lij such that only 



certain configurations of colours are allowed in the statistical sum (2.1) 



In what follows we will be dealing with a lattice formed by coloured faces where 
each adjacent face can not have the same colour. As it was remarked by Baxter in 



this system can be thought as a system of particles interacting through an infinitely 
repulsive force between nearst neighbors of the same type. For the trigonometric SOS 
model we will have = 6 + kj'y, where is an integer variable while 6 and 7 are 
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Figure 5: Face- vertex configurations and the associated Boltzmann weight. 



complex numbers. In the Fig. [5] the height function /j j characterising the colour of the 
face is projected into the center of the j) face. Interestingly enough, it was remarked 
by Lenard an equivalence between this colouring of faces and the configurations of a 
six vertex model. This equivalence is depicted in Fig. |5] where basically one removes 
the outer edges of the four-faces set and places arrows on the internal edges according 
to a certain rule. This rule is as follows: each face of a four-faces set is visited in the 
anticlockwise direction. If the colour changes by +7 when intersecting an edge, this 
edge receives an arrow pointing inwards; and if the colour changes by —7 we place an 
arrow pointing outwards on that edge. For a particular class of statistical weights this 
model is also called Three- colouring model 18 ,20 . In that case the height function /jj 



is conveniently labeled by an element of the ring Z3. More precisely, the Z3 structure 
labeling the colours of the faces is unveiled by considering the remainder after division 
of \j by 3. 



The dynamical Yang-Baxter equation. In 21 23 it was demonstrated that the 



Boltzmann weights of a face model satisfying (2.5) are encoded in the solutions of the 
dynamical Yang-Baxter equation. This dynamical version of the quantum Yang-Baxter 



equation was proposed by Felder in 21 as the quantised form of a modified classical 



Yang-Baxter equation. In that case this modified classical Yang-Baxter equation arises as 



the compatibility condition for the Knizhnik-Zamolodchikov-Bernard equations 24,25 
Previous to that, the dynamical Yang-Baxter equation had appeared in connection to 
the Liouville string field theory in [26] . 

Now let V = f + © f _. be a two-dimensional complex vector space and consider the 
operator TZ{X,6) G End(V(8>V) with X,6 EC The dynamical Yang-Baxter equation for 
the trigonometric SOS model then reads 



7^l2(Al -\2,9- 7/^3)7^l3(Al - A3, ^)7^23(A2 - A3, ^ - 7/^1) 
7^23(A2 - A3, ^)7^l3(Al - A3, ^ - 7/i2)7^l2(Al - A2, 



(2.7) 



where h = diag(l, —1). The Eq. (2.7) is defined in End(Vi 0^2 V3) and the action of 
7^i2(A, 9 — 7/2,3) on vi ^V2^ V3 is understood as 

[n{\,e --fh)vi^V2]^vs , (2.8) 

keeping in mind that h^vs = hv^. In other words, h is simply the eigenvalue of h on the 



"'^See Note added in proof of 19 
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corresponding subspace. The explicit trigonometric solution of (2.7) is given by 



7^(A,( 






\ 





c_(A,e) 






c+(A, 










(2.9) 



with non-null entries 



a±{\d) 



sinh (A + 7) 
sinh (A) 



c±{X,0) = sinh (7) 



sinh (6) 
sinh {9 =F A) 
sinh {6) 



(2.10) 



The solution described by (2.9) and (2.10) consists of a particular limit of the elliptic 



solution found in 21 , 22 . Also it is important to remark here that such solutions are 



in correspondence with Baxter's eight-vertex model after a vertex-face transformation 



27,28j. Moreover, the dynamical 7^-matrix (2.9) satisfies the ice rule 

nabi\,9),ha + h =0 



(2.11) 



which plays an important role in estabilishing an algebra associated to dynamical TZ- 
matrices. 

Now we turn our attention to the relation between solutions of the dynamical Yang- 



Baxter equation (2.7) and the statistical weights of a SOS model satisfying (2.5). At 



this stage the observations of Lenard 19 , the vertex-face transformation introduced 



by Baxter 28 and Felder's dynamical Yang-Baxter equation 22 converge to the same 



point. Following 22 we thus have 

9 9 + ^ 

9 9-^ 

(9--1 9 
9 9-^ 

This association is also depicted in Fig. [5] 



a+{X,9) 




K{\9) 






c+{\9) 





^ + 7 



9 



6/ - 7 

^ + 7 9 + 2-f 
9 9 + -f 

+ 7 9 
9 9 + -f 



a_{\9) 

b-{K9) 
c.{X,9) 



(2.12) 



The dynamical monodromy matrix. Still following 22,23 we define an inhomo- 
geneous monodromy matrix 7^(A, 9) formed by the ordered product of dynamical TZ- 
matrices. More precisely, this dynamical monodromy matrix reads 



ra{\,9)=l[na^{\-^i^,9,) 



(2.13) 



i=l 



7 



with 9i = 9 — -y hk, and it should not be confused with the face monodromy matrix 

k=i+l 

defined in ( [2^ . 

The dynamical monodromy matrix (2.13) is an operator living in the tensor product 
space Vq ® Vi ® • • • ® Vl. The space Ya will be refereed to as auxiliar space while the 
tensor product Vi (8> ■ ■ ■ ® Vl will be called quantum space. In this way Ta{X,9) can 
be regarded as a matrix on the auxiliar space whose entries are matrices living in the 
quantum space. Here we shall restrict ourselves to the monodromy matrix built out of 



the 7^-matrix (2.9). Therefore, it can be conveniently denoted as 



ra{x,9) 



A{X, 
C{X, 



B{X,9) 
D{X,9) 



(2.14) 



Due to the dynamical Yang-Baxter equation ( |2.7 ) and the ice rule (2.11 ), one can demon- 
strate that the monodromy matrix (2.13) satisfies the algebraic relation 



7^a6(Al - A2, ^ - 7H)r„(Ai, 9)%{X2, 9 - jha) = %{X2, 9)Ta{Xr, 9 - ^k)nab{Xi - A2, I 



(2.15) 



where H = hk- With the help of the definition of /i, we can identify H as the Cartan 



k=l 



element of the su(2) algebra on the tensor product space \i ® ■ ■ ■ ® Yl- Furthermore, 
in the limit — t- 00 we can immediately see that (2.9) becomes the standard 7^-matrix 



invariant under the quantum affine algebra ?7g[su(2)]. This observation extends to the 



monodromy matrix (2.13) and the algebra (2.15), which become respectively the standard 



trigonometric six vertex model monodromy matrix and Yang-Baxter algebra 29 , 2] . 



Domain wall boundary conditions. We shall now consider the trigonometric SOS 
model on the lattice described in Fig. [T] with N = M = L + 1 and special bound- 



ary conditions. As for the boundary conditions, we set lij 



I 



^ + (L + 1 - j)7 



and Il+ij = lj,L+i = ^ + (j ~ 1)7- This special boundary conditions is illustrated in 
Fig. |6] together with the corresponding structure of vertices. In the vertices language we 
can immediately recognize this special boundary conditions as the case of domain wall 
boundaries introduced by Korepin in [2] . This observation if of fundamental importance 
allowing us to express the model partition function in terms of the entries of the mon- 
odromy matrix (2.14) analogously to the case of the standard six vertex model |4j. The 
diagrammatical interpretation of (2.13) is given in 22,23 and for a discussion on the 



construction of the partition function (2.1) in terms of the components (2.14) we refer 
to 30 . For the trigonometric SOS model considered here we only have to keep in mind 



that the entries of the dynamical monodromy matrix (2.14) also depends on the colour 
variable 9 governed by the height functions Ij^i. In this way the partition function (2.1) 
for the trigonometric SOS model with domain wall boundaries can be written as 



Ze = {0\l[B{X„9+n)\0) 



(2.16) 
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Figure 6: (a) SOS model with domain wall boundaries. (6) The corresponding 
structure of vertices. 



where 



|0) 



and 



|0) 



(2.17) 



Dynamical Yang-Baxter algebra and functional re- 
lations 



The dynamical Yang-Baxter algebra (2.15) encodes commutation rules for the entries 
of the dynamical monodromy matrix (2.13). In contrast to the standard Yang-Baxter 
algebra, the relation (2.15) not only contains the generators v4(A, 6'), -B(A, 6'), C(A, 9) and 



D{X,6), but also the 5u(2) Cartan generator H. This indicates that the algebra defined 



by (2.15) needs to be complemented. 



The su(2) Cartan generator. The definition of H together with (2.13) and (2.14) 



allow us to directly compute the commutators between the su(2) Cartan generator H 



arbitrary 9. Nevertheless, this analysis can also be found in App. ffland we have 



and the entries of the monodromy matrix (2.14). In the limit 6' — )■ oo this analysis has 
been performed in l2l|4| and here we find that there are no significant modifications for 



(3.1) 



[Ai\,e),R] = 
[C(A,^),H] = -2C(A,i 



[5(A,0),H] = 2fi(A,( 
[D{X,e),E] = 0. 



From the 7^-matrix (2.9) we can see that the generator H will appear in (2.15) only as 



g±7H_ fjj-^^g convenience we define the operator K = with g = e"^, in such a way 



that the commutation rules (3.1) imply the following relations: 



B{X,e)K = q^KB{\,e) 
C{\,e)K = q-^KC{X,e 



[A{\,9),K] = 
[D{\,e),K] = 0. 



(3.2) 



The commutation rules (3.2) in addition to (2.15) form an extended Yang-Baxter 
algebra with generators A{X, 9), B{X, 9), C(A, 9), -D(A, 9) and K^*". For our purposes here 
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we only need to extract a few commutation relations from (2.15). Namely, 
S(Ai,^)5(A2,^ + 7) 



A(Ai,^ + 7)5(A2,^) 



5(A2,0)fi(Ai,0 + 7) 
g(A2- A1 + 7) g(g + 7) 
s(A2 - Ai) s{e + 27) 
s(^ + 7-A2 + Ai) 5(7) 



s(A2-Ai) s(^ + 27) 



S(A2,^ + 7)^(Ai,^ + 27) 
S(Ai,^ + 7)v4(A2,e + 27) 



D(Ai,e-7)5(A2 
s(Ai - A2 + 7) 



5(A2, e - 7)/^(Ai, e)[tqK-' - t-\-^K][tq^K-^ - r^'^K]-' 



s{\i - A2) 
s(Ai - A2 



5(Ai, e - 7)L'(A2, e)[tqxiX^^K-^ - t-\-^x^^X2K\[tq^K-^ - r^q-^K]'^ 



s{9) 



5(A2, e + 7)C7(Ai, 6 + 27) [tqK-^ - t-^g-^K] [tg^R-i - r'q-^K]-^ 
f . ^ • : . ^^ A(A2, e + ^)D{X^, 9)[tqK-^ - t-'q~^K][tq^K-^ - t-^q-^K]-^ 

A(Ai, 6 + 7)£'(A2, e)[tqxiX^^K-^ - t-\-^x:[^X2K][tq^K-^ - t'^-^K]-^ . 



s{9 + ^) 

g(7) 3(^ + 7 + Ai-A2) 
s(^ + 7) s(Ai-A2) 
5(7) 



s(Ai-A2) 



(3.3) 



In the above commutation rules we have employed the notation s(A) = sinh(A), t = 
and Xi = The algebra formed by the relations (3.2) and (3.3) will be one of the main 



ingredients in the derivation of a functional relation determining the partition function 



(2.16). 



In order to proceed we will also need to consider the action of the generators A{X, 9), 



B{X,e), C{X,e), D{X,e) and K on the states |0) and |0) defined in (2.17). Those states 



are the su(2) highest and lowest weight states respectively, and from (2.9), (2.13) and 



(2.14) we readly obtain the relations 



K±|0) = g^^lO) 

L 

A{X,9)\0) = J]s(A-/x, + 7)|0) 



DiX,9)\0) 



(3.4) 



Moreover, we also obtain the properties 

B{x,e)\o) = ] 
B{x,e) |o) = 



C{X,9) |0) = 

c{x,e)\o) = t 



(3.5) 
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where the symbols f and | stand for non-null values. The relations (3.5) together with 



(3.1) support considering B{X,9) and C{X,9) as creation and annihilation operators 



respectively with respect to the pseudo- vacuum state |0). 



Altogether the relations (3.2)-(3.5) allow us to derive the following formula, 



n n— 1 



^^(Ao, ^ + 7) n ^(^- ^ + - 1)^) = Y.^ Il B^rf ' ^ + ^7) |0) 

i=l i=l j=l 

n j—1 n—2 

j=2 1=1 k=l 



(3.6) 



where 



_ J 3 



l<j<i 



j + 1 i ^ j ^ n — 1 



and 



k I < k < i 

k + 1 i < k < j 
k + 2 j <k<n-2 



(3.7) 



In their turn the coefficients Mj and Nji are given by 
M- = 



g(7) s{e + 7) sjXo -Xi + e + {2n-l- L)^) sjO + {n - L)^) sjO + ^7) 
s{Xi - Xo) s{9 + n-f) s{9 + {2n - 1 - L)j) s{9 + {2n - L)j) s{9 + {n - L)-f) 



Y[ s{Xo - + 7)s(A, - III) JJ 



s{Xi - Afc + 7) s{Xk - Aq + 7) 
s(Ai - Afc) s(Afc - Ao) 

t=l fc— 1 

g(7) s(Ao-A, + g + 7) g(g + (n-L)7) g(g + n7) 
■ s(Ao - A,) s{9 + ^7) s(0 + {2n - L)-f) 3(9 + {n - L)-f) 

I. , \ (\ X TT s(Ao - Afc + 7) s(Afc - Ai + 7) 



1=1 



k = l 



s(Ao - Afc) s(Afe - Aj) 



(3.^ 



5(7) 5(7) s(Aj - Xi + 7) s(Ao - A, + ^ + 7) s(Ao - A^ + ^ + (2r2 - 1 - L)7) 



J]^ s{Xi - fj.i + j)s{Xj - Hi) 



s{Xq - Xj) s{Xi - Xq) s{Xj - Xi) s{9 + n-f) 

s{9 + {n-L)-f) s{9 + n-^) 
^ s{9 + {2n - L)7) s{9 + {n - L)^) A 1 

^ TT s{Xj - An, + 7) g(A^ - Aj + 7) 
J-A s{Xj - Xm) s{Xni-Xi) 



s{9 + (2n - 1 - L)7) 
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s(7) 



5(7) s{Xi - Xj + 7) s(Ao -Xj + e + -f) s{Xo - Xi + e + {2n-l- L)-f) 



s(Ao - Xi) s{Xj - Ao) s{Xi - Xj) 



s{9 + ri'y) 



s{e + (2n - 1 - L)7) 



^(W(2.-L)7).(^ + (n-L)7) n ^(^^- - + - 

^ -TT s{Xi - A^ + 7) g(A^ - Aj + 7) 
s(Ai-Am) s(A^-Aj) 



(3.9) 



As expected, the expressions (3.6)-(3.9) recover the ones presented in [4,11 in the hmit 
— ^ 00. At this stage we have gathered most of the ingredients required to obtain a 



functional equation determining the partition function (2.16). 



Functional Equation. The relation (3.6) is valid for the product of an arbitrary 
number of operators B{Xi,6), and following the approach devised in 11 , we examine 
the quantity 



c{Xo,e + ^)l[B{x„d + {z-i)^)\o) 



(3.10) 



i=l 



under the light of the extended Yang-Baxter algebra formed by (2.15) and (3.2). Thus 



considering (3.6) and the definition (2.16) we immediately obtain the relation 

L+l L+1 

(0| C(Ao, ^ + 7) n ^ + - 1)^) |0) = Yl ^i^eiXi, A,_i, A,+i, . . . , A^+i) 

i=l i 
L+l j-l 

+ A'jjZ5i(Ao, . . . , Aj_i, Aj+i, . . . , Aj_i, Aj+i, . . . , Al+i) , (3.11) 

j=2 i=l 



where the coefficients Mj and Nji are the ones given by (3.8) and (3.9) with n = L + 1. 
On the other hand, the highest weight representation theory of the su(2) algebra tells us 



that the LHS of (3.11 ) vanishes. This property has been already discussed in 11 but for 



completeness we also include it in the App. Kj In this way we are left with the following 



functional equation for the partition function (2.16) 

L+l 

M(Ze(Ai, . . . , Aj_i, Aj+i, . . . , Xl+i) 



i=l 



L+l j-l 

NjiZg{Xo, . . . , Ai_i, Aj+i, 

j=2 i=l 



Aj_i, Aj+1, . . . , Xl+i) — . 



(3.12) 



Some comments are in order at this stage. The partition function Zg depends on two 
sets of variables, {Aj} and {fit}, as well as parameters 7 and 6. Nevertheless, within this 
approach we can see that the set of variables {^i} can also be regarded as parameters. 



12 



The Eq. (3.12) is linear and homogeneous in Zq and this observation will have important 
consequences for the characterisation of the desired solution. The homogeneity of (3.12) 
tells us that if Zq is a solution, so is uZq where a is a constant. In fact a only needs 
to be independent of the variables {A^}. Therefore, the Eq. (3.12) can determine the 
partition function at most up to a constant and the full determination of the partition 
function will require that we are at least able tocompute Zq for some particular value 
of the variables Aj. Moreover, the linearity of (3.12) tells us that if Zr^^ and zf^ are 



two solutions of (|3.12l), the linear combination Zf'' + zf^ is also a solution. 



This fact 



suggests that classifying the classes of unique solutions of (3.12) is an important step in 
this framework. The Eq. (3.12) has the same structure of the functional equation derived 
in |Tl], however the coefficients Mj and Nji are deformed by the dynamical parameter Q. 
In what follows these issues will be discussed and the desired solution of (3.12) will be 
presented. 



4 The partition function 

In this section we shall consider extra properties expected for the partition function 



(2.16) which will enable us to select the appropriate solution of (3.12). These properties 
include the multivariate polynomial structure of the partition function Zq^ as well as its 
asymptotic behaviour. In order to avoid an overcrowded section we shall discuss the 
mentioned properties in the App. |B]and here we present only the required results. 

Polynomial structure. The partition function Zq is of the form, 

Z,(A,.....A,) = %^^. (4,1) 



i=l 

where Zq[x\^ . . . ,xl) is a polynomial of order L in each variable Xj = e^'*'' separately. 
The same polynomial structure holds if we consider the variables {/ij} instead of {Aj}, 
although this property will not be required. 

Asymptotic behaviour. In the limit Xi — ?■ oo, the function Zg possesses the asymp- 
totic behaviour 

Z {xi...xl) , (4.2) 

n=l 

where = 1(1 + g)(l + g + g^) . . . (1 + g + ■ ■ ■ + g") denotes the g-factorial function 

and Uj = e^^\ 



For the moment we leave the properties (4.1) and (4.2) at rest and proceed with a 



more careful examination of the functional equation (3.12 ). Firstly we notice that besides 



the set of variables {Ai, A2, . . . , Xl] required to characterise its solution, the Eq. (3.12) 
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also depends on variables Aq and A^+i. Thus the variables Aq and A^+i can be fixed in 
order to fulfill our needs. In particular, those variables can be chosen in such a way that 
solving (3.12) under the conditions (4.1) and (4.2) becomes systematic and simple. For 
illustrative purposes, let us see how this approach would work in the case L = 2. In that 
case our functional equation reads 



MiZe{\2, As) + M2Zg{\,, A3) + MsZg{\,, A2) 

+N2iZg{Xo, A3) + NsiZeiXo, X2) + A^32^e(Ao, Ai) = 



and we set Aq = /ii and A3 = /ii — 7. By doing so we find 



Ml 



->>0=M1 



M, 



2 I ^0 = A'l 
^3=^1-7 







(4.3) 



(4.4) 



and we also define 



m2= M3I Ao=Mi 

A3=A»i-7 
A3=/il-7 



m2 = Ns2\ Ao=Mi 

A3=/il-7 



s{9 + -f) 



sie + 37) 

s(^ + 7 - Ai + /ii) 



S{lfs{fli ~ fl2+ 7)s(/i2 - /il + 7) 



X 



s{e + 37) 

s(A2-/ii) s(A2-Ai + 7) 



s(7)^s(/i2 - /xi + 7)s(Ai - /i2 + 7) 



s(A2-/ii+7) s(A2 - Ai) 

s(e + 7-A2 + /ii) , -2 / , \ f\ , \ 
^ s(7) s(/i2 - /ii + 7)s(A2 - /i2 + 7) 

g(Ai-/ii) g(Ai-A2 + 7) 
s(Ai-;Ui+7) s(Ai-A2) 



X 



(4.5) 



With the above specialisation of the variables Aq and A3, the Eq. (4.3) reduces to 



, . s(^ + 7 - A2 + /ii) s(A2 - /i2 + 7) s{\i - /ii) s(Ai - A2 + 7) 7 . , ^ 

^e(Al,A2) = 7^- r ^ r-TT ■ r tt ^-^ Zg{fli, Xi) 



+ 



5(6* + 7) s(/ii - /i2 + 7) s(Ai - /ii + 7) s(Ai - A2 
s(6' + 7-Ai+/ii)g(Ai-/i2 + 7) g(A2-/ii) g(A2-Ai+7) 
s{9 + 'j) s(/ii - /i2 + 7) s(A2 - /ii +7) s(A2-Ai) 



2'e(/ii, A2) 
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M3 



A3=Ml-7 



Zg{^l,fil - 7) 



(4.6) 



In addition to that the Eq. ( |4.6[ ) reduces to the following identity when A2 = /^i, 

"iV2l" 



•''0'-"'2=Ml 
A3=A'i-7 



Zg{ij,i, /ii - 7) = Z0(/ii, Ai) - Ze(Ai, /ii) 



(4.7) 



As demonstrated in the App. |D| the functional equation (|3.12|) admits only symmetric 
solutions, i.e 

the quantity 



solutions, i.e. Zg{Xi,X2) = Zg{X2,Xi), and consequently the RHS of (4.7) vanishes. As 

'N2 



M3 



Ao,A2=Mi is finite we can thus conclude that Zg^iii.^i — 7) = 0. This 

-^3=^1-7 
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property simplifies (4.6) and we are left with 



- fi2 + l)s{0 + 'y)Zg{Xi, A2) 

s(6' + 7 - A2 + /xi)s(A2 - /i2 + 7) 



s(6' + 7 - Ai + /ii)s(Ai - /i2 + 7) 



g(Ai-/ii) g(Ai-A2 + 7) 
s(Ai-/ii + 7) s(Ai-A2) 

g(A2-/ii) g(A2-Ai+7) 
s{\2 - fJ'i + 1) s(A2-Ai) 



Ai) + 
^e(/^i, A2) . 



(4i 



The vanishing condition of Zq above unveiled have a special appeal since we are 
interested in the polynomial solution of (4.3) with order dictated by (4.1). For instance, 
they allow us to write 

ZeiXi, X2)\,,=,, = s(A2 - /ii + 7)^(A2) (4.9) 

where now V{X2) needs to be a polynomial of the same order as Ze with L = 1 in order 
to satisfy (4.1). The expression (4.9) can now be replaced in (4.8) yielding the following 
relation, 



Ze{Xi, A2|/ii, /i2) — 

s(6' + 7 - A2 + jUi) ■s(A2 - /i2 + 7) 

s(6' + 7) s(/ii-/i2 + 7) 
5(6* + 7 - Ai + fii) g(Ai - H2 + 1) 

s(6' + 7) s{ni-H2 + 7) 



M ^ 3(Ai-A2 + 7) ,^^^ ^ , 

s(Ai - /ii) — T- — — V{Xi) + 



s{X2 - /il 



six, - A2) 
g(A2 - Ai + 7) 
s(A2-Ai) 



V{X2 



(4.10) 



which can be substituted back into the original equation (4.3). This step leave us with 
an equation involving the functions V{Xo), ^(Ai), l^(A2) and ^^(As) for arbitrary values 
of those variables. Moreover, by setting A3 = /ii we are left with the equation 

K,V{X2) + i^2^(Ai) + L2iV{Xo) = , (4.11) 

where the coefficients Ki, K2 and L21 coincide respectively with Mi, M2 and N21 obtained 
from (3.8) and (3.9) with L = 1 {n = 2), 6 ^ 9 -\- j and fii — > fi2- So the function V{X) 
obeys the same equation as Zg+j{X\fi2)- 

We have now reached an important stage of this approach. Let us suppose that the 
solution of (3.12) with polynomial structure (4.1) is unique. In fact, the uniqueness 



of the polynomial solutions is demonstrated in the App. [Ej Since V^(A) and ^^(A) are 
polynomials of the same order, that is to say V"(A) = f]2Zg+^{X) where i72 does not 
depend on A. The solution of (3.12) for L = 1 can be found in the App. [F] and here 
we shall only make use of the solution. In this way the results so far can be gathered 
and from (4.10) and (F.5) we immediately obtain an explicit solution for Zg{Xi, A2). The 

Ze{Xi, A2) = F12 + F-^ 



solution is then given by 
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(4.12) 
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where 



s{9 + - Xi + fii) s{9 + 2-f - Xj + fi2) s{Xi - H2 + l) s{Xj-Xi + -f) 



+ 7) 



s{e + 2-f) s(/ii-/i2 + 7) 



s{Xj - Xi) 



(4.13) 

The constant factor can be fixed by the asymptotic behaviour (4.2) and we find 



= ■s(7)^s(/ii — /i2 + 7)- Thus we have completely determined the partition function 



2.16) for L = 2 using solely the polynomial structure (4.1) and the asymptotic behaviour 



(4.2), in addition to the functional equation (3.12). In what follows we shall consider the 



general L case. 

Special Zeroes and Symmetry. The first step in order to consider the case with 



arbitrary values of L is to obtain an analogous of the relation (4.9 ) which can be obtained 



by uncovering special zeroes of our partition function. These zeroes have been unveiled 
in App. O and in addition to that we shall also make use of the symmetry property 



discussed in App. |D| Thus taking into account (C.13) and (D.5) we can write 

L 

i=2 



(4.14) 



We proceed by setting Aq = /ii and A^+i = /ii — 7 in the functional equation (3.12). By 
doing so we obtain the expression 

L L 



Ze{Xi, . . . , Al) = TT s(Afe - Ail + 7)^\/(Ai, . . . , Aj_i, A^+i, . . . , A^) 
1=1 fc=l ^ 



where 



(4.15) 



rriL 



rrii 



s{9 + j) 



s(^ + 7 - Xj + /ii 
s(0 + (L + l)7) 



'l)^^'s(7)2 II s(/ii - fij + 7)s(/i, - /ii + 7) 



J=2 
L 



-1)^5(7)^ n ^(^'^ - /^i + 7)s(Aj - /ifc + 7) 



k=2 



X 



g(Afc-/ii) g(Afc-Aj+7) 

s{Xk - /ii + 7) s{Xk - Xj) 
+3 



n 



(4.16) 



It is important to remark here that we have also considered (4.14) and the symmetry 



property Zq{. . . , Xi, . . . , Xj, . . .) = Zg{. . . , A^, . . . , Aj, . . . ) discussed in the App.[D|in order 
to obtain (4.15). 
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The relation (4.15) can now be substituted back into the Eq. (3.12) and considering 
Xl+1 = A^i we obtain 



KiV{Xi, . . . , Aj_i, Aj+i, . . . , Xl) 



i=l 



L j-1 



+ LjiV{Xo, . . . , Ai_i, Aj+i, . . . , Aj_i, Aj+i, . . . , Xl) — . 

j=2 i=l 



(4.17) 



In their turn the coefficients Ki and L^j appearing in (4.17) correspond respectively to 
the coefficients Mj and Nji given in (3.8) and (3.9) with L ^ L — 1 {n = L), 6^9 + '-/ 
and /ij — )■ /ij+i. Moreover, the compatibility between (4.14) and (4.1) tells us that the 
function is a multivariate polynomial of the same order as the partition function Zg 
for a lattice with dimensions {L — 1) x (L — 1). Thus the Eq. (4.17), together with the 
uniqueness property discussed in the App. |E| implies in 

V{Xi, ...,Xn) = ^nZe+^iXi, . . . , A„) . (4.18) 

In this way the relation (4.15) can be iterated using the results obtained in the App. [F] 
as initial condition. 

By carrying on with this procedure we obtain the following expression for our partition 
function: 

Ze(Ai,...,Ai)= Yl ^^1-^ (4.19) 

where 



ii...tl 



Uk=2 -f^k + l) „=i 



+ ?27 - Ai„ + /in) -A- -A- \TT r\ 
s(e + n^) ii ii ^(A.„ -/i,+7)lls(A.„-/x, 



n=l j=n+l 



X 



L-1 L 

nn 

n=l m>n 



g(Ai^ - +7) 



(4.20) 



Here Sl denotes the permutation group of order L and the asymptotic behaviour (4.2) 
implies in i^^ = s{'y)^ Ilt=2 ^if^^ - /^fc + t)- 



5 Multiple integral representation 



The expression ( |4.19[ 4.20) can be converted into a multiple contour integral similarly to 
the expression recently found in 30 for the ?7q[su(2)] vertex model. As a matter of fact. 



multiple contour integrals seems to fit naturally into the algebraic-functional framework 



presented here. We start by noticing that the solution of (3.12) for L 
can be rewritten as 



1 given in (F.5) 



Ze(X) 



27ri / s{w — A) 



+ 7- w + /ii) 
s(0 + 7) 



dw 



(5.1) 
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where the integration contour contains the pole at w = X. Now we look to the Eq. (4.15) 



considering (4.18) and keeping in mind that for L = 1 we have (5.1). This suggests that 



the iteration procedure described by (4.15 ) can be mimicked by Cauchy's residue formula. 



It turns out that when we look to (4.15) searching for solutions as contour integrals, we 



find a factorised formula for the integrand. In this way we end up with the following 
expression for our partition function, 



Zg{Xi, . . . , Xl) 



s(7) 
27ri 



nf=in 



{wj -Wi + 'y)s{wj - Wi) 



n 



+ jl - Wj + fij) 



Ui,j=lS{w^- Xj) 
L i-1 L 

n n ^^^^ " n ~ + ^^^^ • • • ^^^^ 

i=\ j=l j=i+l 



X 



(5.2) 



where the integration countours enclose the poles at Wi = Xj. As expected the expression 



(5.2) coincides with (4.19)-(4.20) when evaluated using Cauchy's residue formula. More- 



over, in the limit 6 ^ oo the formula (5.2) reduces to the one obtained in |30] after a 
relabelling of the variables fij. This relabelling does not affect the solution Zq{X\^ . . . , Xl) 
since this partition function is invariant under the exchange of variables fj,i -H- fj,j as dis- 
cusssed in [ij. 



6 Concluding remarks 



The main result of this paper is the integral representation (5.2 ) obtained for the partition 
function of the trigonometric SOS model with domain wall boundaries. This integral 
formula has been obtained by solving a functional equation derived from the dynamical 
Yang-Baxter algebra. This approach has been proposed in 
present a more robust formulation of that method. 



11,12 and here we also 



In contrast to the case considered in pT,12 , where the su(2) algebra only appears in 
the final stages of the derivation of (3.12), here it plays an important role from the very 



(3.2) and (3.4) 



beginning. For instance, the derivation of (3.6) requires the repeated use of the relations 



It is important to remark here that the elliptic version of this same partition function 
has been considered previously in l7]-[9,31 



In particular, the work [8] discusses the lack 
of a single determinant expression for this partition function generalising the Izergin- 
Korepin determinant. For the three-colouring model functional equation for 



this partition function was obtained in 32 , 33 though a connection with the functional 



equation presented here is not obvious at the moment. It is also worth remarking that 
the trigonometric SOS model with one reflecting end, and the remaining boundaries of 
domain wall type has been considered in 34 ,35 . In that case the dynamical Yang-Baxter 



algebra also plays an important role, though it is only responsible for a few out of six 
conditions determining uniquely the model partition function. The approach considered 
here makes use of only three conditions and it would be interesting to extend it to the 
case considered in 34 , 35 
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Moreover, it has been recently discussed in |36j the usefulness of such integral formulas 
for computing correlation functions for the case of domain wall boundaries which makes 
the representation (5.2) more attractive. The generalisation of our results for the elliptic 
case is under investigation and we hope to report on that in a future publication. 
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A Dynamical Yang-Baxter algebra vs. su(2) 



The analysis performed here will follow the same lines as the one presented in [12]. We 
shall consider the 5u(2) generators E, F and H satisfying the relations 



[E,F]=H 



[H, E] = 2E 



[H, F] = -2F 



(A.l) 



and compute their commutation relations with the generators of the dynamical Yang- 
Baxter algebra A{X,9), B{\,9), C{X,9) and D{X,9). In fact we will only need their 
commutation rules with the Cartan generator H whose fundamental representation on 
the quantum space is given by 



(A.2) 



1=1 



Here hj consists of the Pauli matrix 



h 



1 
-1 



(A.3) 



acting non-trivia lly on the space Vj of the tensor product Vi ® ■ ■ ■ ® V^. 

The ice rule ( 2.1l[ ) can be rewritten as [Tlaj{X,9j),hj] = —[JZaj{\,9j),ha] which im- 
mediately lead us to the relation 



[Ta{X,9),}i] = -[ra{X,9),k 



(A.4) 



due to the definition (2.2 ). In terms of the monodromy matrix entries (2.14), the relation 
(A.4) exphcitly reads 



[A{X,9),E] = 
[C{X,9),11] = -2C{X,^ 



[B{X,9),E] = 25(A,( 
[D{X,9),R] = 0, 



(A.5) 
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which allows us to exploit the representation theory of the su(2) algebra in order to gain 
insight into the dynamical Yang-Baxter algebra generators. 

For instance, the su(2) highest and lowest weight states |0) and |0) defined in (2.17) 
obey the relations H |0) = |0) and H |0) = — L|0). These properties together with 
(A. 5) allow us to obtain the relation 

n n 

H J] B{K, e + {i- 1)7) |0) = (L - 2n) J] 5(A„ 9 + ( 



1)7) |0) 



(A.6) 



which is valid for any number n of operators B{X,6). Now the relation (A.6) put us in 
position to use the su(2) representation the ory t o draw conclusions about the generator 
B{\, 9). For the case n = L the expression (A.6) tells us that Y[i=i 9+ {i — 1)7) |0) 
is an eigenvector of H with eigenvalue —L. On the other hand this is the same eigenvalue 
associated with the state |0). Since this eigenvalue is not degenerated we can conclude 
that 



l[B{\,9 + {z-lh)\0)^\0) 



and from (3.5) we immediately have that 

L+l 

ni?(A^,^^ + (^-l)7)|0) =0 



(A.7) 



The property (A. 8) is an important ingredient for the derivation of the functional equa- 
tion ([3T2I). 



B Polynomial structure and asymptotic behaviour 



In order to analyse the dependence of Zg with the set of variables {Aj} we first consider 
the following change of variables: 



Xi 6 



Xi 6 



2A, 



Ui = e 



q = e' 
t = e' 



(B.i: 



In terms of the above defined variables, the 7^-matrix given by (2.9) and (2.10) can be 
written as 

-(xU + V) , (B.2) 



7^ 



where 
U = 



+ 



(3g' + l)l®l + (g'-l)H®H + 

^ VE O F + —r^ r^F O E 



qx 

(g2-l)(t2 + i) 



it' - 1) 



[t^ - 1) 



it' 



it' 



V 



■ (3 H- 8 IH- - 1)H 8 H - J^''^ ^ 1 « H H- J^''^ « 1 

it' - 1) it' - 1) 



+ 



At'(q'-l)^ ^ 4(1 



it' 



it' - 1) 



(B.3) 
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In (B.3) the generators E, F and H are the su(2) generators satisfying (A.l), and con- 
sidering (B.2), (2.13) and (2.14) we readly obtain the expansion 



B{Ke) 



(B.4) 



Now looking to the product ni=i ^i^j^ ^ + J7) appearing in the definition (2.16), we 



can conclude that 



Zg{\i, . . . , Xl) 



Ze{xi, 



.XL 



(B.5) 



n 



X, 



where Zg{xi, . . . , xl) is a polynomial of order L in each variable Xi. 

Also from (B.2 ) we can see that in the limit x — )■ oo only the operator U contributes for 
the partition function Zq. In (B.3) the operator U is written in terms of su(2) generators 

ins we 

(B.6) 



which allows us to follow the same analysis of 11 . Without significant modifications we 
find that 



2L2 



-{xi . . .Xl) as Xj — )■ cxD . 



n=l 

Here the g-factorial function is defined as 

[n],\ = 1(1 + g)(l + g + g2) . . . (1 + g + . . . + g«) 



(B.7) 



C Special Zeroes 



One important ingredient for solving the Eq. (3.12) under the conditions (4.1) and 



(4.2) is the localisation of some special zeroes of our partition function. Since we are 



interested in the polynomial solution of (3.12), those zeroes will play an important role 



in the characterisation of our solution. We shall start by looking to particular values of 
L for illustrative purposes and then we proceed to the general case. 

• L = 2: 

We set A3 = /ii and A2 = /ii — 7 in such a way that the functions M2, M3, A^2i, -^31 and 
A^32 vanish. For these particular values of A3 and A2 we are thus left with 



Since Mi\^^ is finite we can conclude that Zg{fii — 7,/ii) = . 
• L = 3: 



{c.i: 
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By setting A4 = /xi and A3 = /ii — 7 we obtain 



i=0 



where 



21IA3,A4 



, Pi = M2|;,3 and P2 = Ml 



(C.2) 



(C.3) 



In terms of the variables Xi, the functions Pi are rational functions and thus 3 Aj : Pj = 0. 
Besides the above specialisation of the variables A4 and A3, we also choose A, | Pj = 
for i = 1,2. Thus we are left with 

Po|,^,^Z,(/ii,/ii-7,Ao) = 0, (C.4) 

and since Po\x^ ^2 fiiiite we can conclude that Zg{fii,fii — 7, Aq) = 0. 

• L = 4: 

For the case L > 4 this analysis becomes a bit more involved. We start by setting 
A5 = fii and A4 = ;Ui — 7 similarly to the previous cases. Under this specialisation the 



Eq. (3.12) reduces to 



^llA4,A5^e(A2,A3,/il -7,/^l) + ^2!;,^ ;^^ Z0(Ai, A3,/il -7,/^i) + 
^Zg{Xi,X2,Hl--f,Hl) + N2l\. ^ Ze(Ao,A3,yUi -7,/ii) + 



3IA4,A5 



^3iIa4,A5 ^^(-^O, A2,/il -1,^1) + iV32|A4,A5 ^e(Ao, Ai,/ii -7,/^i) = 



(C.5) 



Next we set Aq = fJ-i and Ai = /Xi — 7. The Eq. (C.5) does not suffer significant 



simplifications and we then proceed by setting A3 = /ii using the following properties: 

Ml 

Ao,Ai,A4,A5 

= -1. (C.6) 



lim 

Aa-^/ti A'31 
M2 



lim 



By doing so we end up with the relation 



Ao,Ai,A4,A5 



M 



3lAo,Ai,A3,A4,A5^e(/il-7,A2,/il-7,/il) = N^. 



21 1 Ao,Ai,A3,A4,A5 



2'e(/ii,/ii,/ii -7,/ii) , (C.7) 



which is further simplified to 

^2l|Ao,Ai,A2,A3,A4,A5^e(/il,/il,/il -7,/^l) = (C.8) 

is finite which allow us to conclude 



with A2 = /ii - 7- The quantity A^aij^^^^^^^^ 
that 

Ze(/ii,/ii,/ii - 7,/ii) = . 



(C.9) 
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Now we move backwards considering the consequences of (C.9) to the previous equa- 
tions. From (C.8) and (C.9) we have that 



(C.IO) 



Ai. This yields the following expression 



which can be reintroduced in (|C.5|) with the above mentioned specialisation of Aq and 



Ze(A2,A3,/ii -7,Aii) 
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Ml 



Ml 



Ao,Ai,A4,A5 



Ao,Ai,A4,A5 



Z0(/ii,A2,/ii -7,Ati) • (C.ll) 



Now we replace (C.ll) back into (C.5) to obtain the expression 

3 



^ Qi Ze{jji, Ai, /ii - 7, /ii) = . 



(C.12) 



i=0 



The explicit form of Qi is not enlightening and shall not be presented here. Nevertheless, 
using similar arguements as for the cases L = 2, 3 we can conclude that Ze(/ii,A,/ii — 



7,/ii) = 0. Thus from (C.ll) we obtain the vanishing condition 

Ze{X2, A3,/ii - 7,/ii) = . 

• General L: 



(C.13) 



For arbitrary values of L we initially set A ^+1 = /ii and A^ = ;Ui — 7 in the functional 
equation (3.12), followed by the specialisation Aq = /xi and Ai = /xi — 7. We collect 
the results at each one of the steps and then start fixing the variables Xl-i = /^i, 
Al-2 = fJ'i — 1, Al_3 = yUi and so on until we exhaust all the variables. Then the 
consistency condition of each step with the previous ones allow us to conclude that 



Zg{fii,fli - 7, A3, ... , Al) 



(C.14) 



Together with the symmetry property discussed in App. D, the relation (C.14) plays an 



important role for solving (3.12) 



D ^0 as a symmetric function 



In this appendix we intend to show that Eq. (3.12) admits only analytic solutions which 
are symmetric under the exchange of variables Aj O Aj. This is an expected property of 



our partition function (2.16) due to the commutation relations (3.3). Nevertheless, we 



shall demonstrate that this property is not an extra input required to solve Eq. (3.12). 
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We start by integrating the Eq. (3.12) over the contour Cj containing only the 
variable Xj. For a given j, the coefficients Mj and N^i {k,l ^ j) do not contain poles 
when Ao — > Xj- Moreover, we also have the following identities between the coefficients 



lim s(Ao — \j) Mk 
lim s(Ao — \j) Mk 

Xq-^Xj 



lim s(Ao - Xj) Njk 
lim s{Xo- Xj) Nkj 

Aq— >A, 



k<J 
k>J 



for j = 1,. . . ,L. Thus after the integration of (3.12) over the contour Cj, we are left 
with the relation 



E 

i=l 
^3 



M, 



L+l 



[ZeiX,,.. 

Zg{Xj, Ai, 



where 



Aj_i, Aj+i, . . . , Al+i) — Zg{Xj, Ai, . . . , Ai_i, Aj+i, . . . , A^+i)] — 
..,XL)-Ze{X,,...,XL) (D.2) 



Mfc = lim s(Ao - Xj) Mk . (D.3) 

Aq— >A, 



Two observations are important at this stage. Firstly, the relation (D.2) is valid for 
j = 1, . . . ,L and thus it provides us with a total of L equations. Secondly, we notice 
that the RHS of (D.2) does not depend on the variable Al+i. In fact this variable can 
be adjusted, together with the results obtained for the (j — l)-th equation, in order to 

(D.4) 



show that the RHS of (D.2) vanishes. Thus the relation (D.2) implies in 



Zg{Xi, . . . , Xl) — Zg{Xj, Ai, . . . , Aj_i, Aj+1, . . . , Xl) 

for j = 1, . . . , L and consequently we have the desired symmetry relation 
Zg{Xi, . . . , Xi, . . . , Xj, . . . , Xl) = Zg[Xi, . . . , Xj, . . . , Xi, . . . , Xl) 



(D.5) 



E Uniqueness 

In this appendix we prove the uniqueness of the multivariate polynomial solution of the 
Eq. (3.12). In order to start we ffist need to introduce the modified coefficients 

L+l L+l 

m = Y[ xJ^Mi and Nji = Y[ x'^Nji , (E.l) 







with Xi = e \ In this way the Eq. (3.12) is given by 

L+l 

MiZg{Xi, . . . , Aj_i, Ai+1, . . . , Al+i) 



i=l 



L+l j-1 

+ NjiZg{Xo, . . . , Aj_i, Aj+1, . . . , Aj_i, Aj+1, . . . , Ai+i) = 

j=2 i=l 



(E.2) 
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in terms of Zq which is a polynomial of order L in each variable Xi according to (4.1) 



We shall now explore the linearity of the Eq. (E.2). More precisely that means the 



following: if Z\ and Z2 are two multivariate polynomials of type (4.1) satisfying (E.2), 
then 

Z = aZx- [iZ2 (E.3) 

is also a solution for any constants a and (3. Polynomials are characterised by the location 
of their zeroes and we can express Z, Zi and Z2 as 



X — Tj 



X Si 



i=l 



Zo 



11(^-^0 

i=l 



(E.4) 



where x can represent any of the variables x,. Since the constants a and (3 in (E.3) are 



arbitrary they can always be fine tunned in order to ensure that Z is also of order L. 



Next we set x = rj in (E.3) and from (E.4) we obtain 



a 



1=1 



Tj Si 



1=1 



(E.5) 



The relation (E.5) allows us to make important conclusions. For instance, if we assume 
that {rj} 7^ {si} then (E.5) implies that {sj} = {ti} since a and (3 can always be 



adjusted to compensate an overall factor. This implies that Zi and Z2 are proportional 



to each other and so is Z due to (E.3). This consequence clearly contradicts the initial 



assumption {vi} 7^ {si}- The remaining option is allowing {rj} = {sj} and thus Z 
and Zi only differ by a constant. This fact together with (E.3) tell us that Z2 is also 



proportional to Zi. In summary this analysis shows that if we have two polynomials of 



the same order solving (E.2), they are essentially the same polynomial. This proves the 



uniqueness of the polynomial solution of (3.12) 



F Solution for L = 1 



Here we shall present the solution of the Eq. (3.12) for the case L = 1 which is of 
fundamental importance in order to derive the solution for general L. For the case L = 1 
the Eq. ([3l2| reads 



MiZg{X2) + M2Ze{\i) + N2iZe{\) = . 



(F.i: 



At first look the condition Aj = A does not seem helpful in finding the solution of (F.l) 



However, a closer look reveals that the coefficients Mi, M2 and A'"2i contain poles when 
the variables Aj coincide. As we shall see this fact will be of fundamental importance. 
We can compute the limit Aj = A of the Eq. (F.l ) using L'Hopital's rule and we are 



left with the following second order differential equation 



PoZe + Pi 



dZe 
dx 



P2 



d^Ze 
dx^ 



(F.2) 
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given in terms of variables x = e^^ and u = e^^^. The coefficients in (F.2) are given by 



Po = (-4g^ + 2qH^ + 2qH')u + (2g^ + 2g* - AqH')x 

P, = {-4qH^ + 2qH^ + 2qH^)u^ + (4^^ - Aqh^)xu + (-2^^ - 2q^ + 4g6t')x2 



AqH' + qH^ + qH*)xu 



and by standard methods we find the general solution 

Zq{x) =Ci{x~ qH'^u) + C2{x- qH'^u) 



_4g2 _ q2^2 ^ 5^4^2 ^ 5^6^2 _ ^8^2 _ 4^8^4)^^2 

(F.3) 



(x — qH'^u) 



dx , 



(F.4) 



where Ci and C2 are two arbitrary integration constant s. N ow the polynomial structure 

(g-g-^)n 



(F.5) 



(4.1) asks for C2 = 0, while the asymptotic behaviour (4.2^ implies in C 
— J — 

t'^q'^)^^Ui ^. Thus our partition function for L = 1 is given by 



+ 7 - A + /ii) 
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